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Abstract 

The 2d one component gas of pointlike charges in a uniform neutralizing background interacting 

with a logarithmic potential is a common model for plasmas. In its classical equilibrium statistics 

at fixed temperature (canonical ensemble) it is formally related to certain types of random matrices 

with Gaussian distribution and complex eigenvalues. In this paper, I present an exact integration of 

this ensemble for N such particles (or alternatively N x N matrices) for all complex temperatures, 

a significant open problem in statistical physics for several decades. I achieve this exact integration 

■ via an exact integration of a related ensemble, the two-dimensional Selberg integral. 

Q-f 
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Non-trivial exact solutions to statistical ensembles with finite degrees of freedom are rare. 
The complexity of these ensembles normally requires perturbation or mean field approaches 
or Monte Carlo simulations to extract meaningful results. Examples of fixed temperature 
(T = 1/ ft) ensembles with exact solutions are the two-dimensional two-component Coulomb 



gas 



2(j and t 



;he one- dimensional one-component Coulomb gas with harmonic confinment 
potential [7j][6] [l6l. The two-dimensional one-component Coulomb gas has been solved 
exactly for f3 = 2 [s| and up to N = 4 jl3| with speculation that it may be exactly solvable 
for even integer temperatures J2l| but no definite result. Pair correlations have also been 
studied numerically at f3 — 4 [22| and compactness has been studied in [15J. The N- vortex 



problem, which is mathematically related given conservation of angular momentum 3j, is 
also considered in fis| . 

The two-dimensional, one-component plasma (OCP) model can be thought of as system 
of N point charges, vortex filaments or lines of charge in the complex plane interacting with 
a logarithmic potential. Each charge has the same sign and value, scaled to be 1, and a 
neutralizing background with charge N. An axisymmetric harmonic potential confines the 
charges towards the origin, giving the energy functional, 

if = ^N 2 -^io g | Zl - Zj |, (i) 

i=l i<j 

where zi is the position of charge % in the complex plane. The statistical mechanics of these 
charges are given by the probability density, 

p(z 1 ,...,z N ) = (Z N ((3))- 1 e^ w , (2) 

where the normalizing constant, 

Z N {p)= [ z 1 zl---z N z* N e-? w 1 (3) 



n 

is the partition function [2j. This particular statistical ensemble was first proposed in 1965 
and solved at the time for one special case of /3 = 2 jg] and solutions for other f3 have 
remained open problems. 

Formally, this statistical ensemble corresponds to the distribution of eigenvalues of several 
classes of complex random matrices where the elements have Gaussian distributions [li]] js|. 
Random matrices have important applications to the quantum theory of atomic nuclei. 
Energy levels are the eigenvalues of Hermitian or unitary random matrices with large order 



N. Because the matrices are Hermitian, the eigenvalues are restricted to the real line. In 



other applications such as vortex lino pinning in superconductors lp| 



121 ] and growth 



models in population biology [17|, however, random matrices are often non- Hermitian and 
may have complex eigenvalues. Feinberg and Zee have made a study of non-Hermitian 



random matrices (j]|5| related to the vortex pinning problem. Ot 



rer applications of non- 



iermitian random matrices are QCD at finite chemical potential 
l^jlfl- For a review of random matrices in quantum physics see (9]. 



19| and quantum chaos 



In this paper, I give the exact, closed form expression for for all even integer (3, 

Z K = W " r -- M f= a F lv (/?),, (4) 



where for = 2 A we have, 
F N (2X) = 



n e 

a<k<l<N pti=0; 



r(l + A) 



r(i+ Pij )T(i + qij ) 



Cp,NCq,N (i 



(5) 



where pki + qui = A and (Yli<k<i<N J2p kl =oj indicates iterated sums rather than a product 
of sums and we have the following recurrence relation: 



n-l 

-r,n (fl j)Cr,n— j — 1 J J 
k=j+l 



r(2(r wfc + YU rja) + \Ea< b < k 2ra & Al) 
r(\E a<b < k+1 2r ab /(k + l)]) 



(6) 



if r nl , 



, i nj 



and rki — for all k — 2, . . . , j + 1 and I — 1, . . . , j, and j = if r n \ ^ 0. 



The [•] means the smallest integer greater than •. This recurrence I will not write out 
in full formula because of its complexity and because this form is most directly related to 
computing it. It is a closed form computation polynomial in N. Note that c rj „ = 1 for 
n < 1. 

I begin by looking at a related integral coined the two-dimensional Selberg integral (not 
to be confused with the complex Selberg): 

TV 



5jy(Ai, A2, A3, A4, A) 



'[o,r]2Jv 

j I [Ofc - xi? + (Vk ~ Vlf] 



n 



x 



(7) 



KkKiV 
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By the change of variables Xj — > (1 — Xj/L)/2 and yj — > (1 — yj/L)/2, we have the partition 
function from the 2D Selberg integral, 

Z N (2\) = f2X)-^'"C-^ x 

li,n4-( 2i )~.^(f4,f,f,A), (8) 

which is a standard way of converting Beta into Gaussian distributions. 
The two-dimensional Selberg integral is a sum of Beta functions. Since 

r(i + Ai)r(i + A 2 )r(i + A 3 )r(i + a 4 ) 

T(2 + Ai + A 2 )r(2 + A 3 + a 4 ) 
= f xjx x ^(l-x) X2 (l-y) x \ (9) 



and 



Kk<KN 



N 



^ 1 c ni,mi,...,njv,mjv J_ J_ j V. 

l<m,mi,...,i 

we have the evaluation, 



l<n±,mi,...,n N ,m N <N(N-l)\ j=l 



Cni,mi,...,rijv,mjv ^ 



sn(^i, a 2 , A3, A4, a) — 

l<ni ,mi,...,npf,mff<N(N— 1)A 

r(i + Ax + 7ij)r(i + A 2 )r(i + a 3 + m J )r(i + a 4 ) 

AA r(2 + Ai + A 2 + nj)Y(2 + A 3 + A 4 + mj) ' 1 j 

where the coefficients c ni>mi ^ nN ,m N are unknown. 

Lemma 1. The values of coefficients c ni)nili ... i „ JViTOjv are independent of order of 
(ni, mi), . . . , (njy, mj\r). -For < 77,1 + mi < n 2 + m 2 < • ■ ■ < 77^ + mjv, the non-zero 
Cni,mi,...,nisr,m N occur when ni,mi, . . . ,n N ,m N satisfy Y^j=i n i + m i = N(N — 1)A and i/ie 
inequalities (j — 1)A < Uj + nij and nj + mj < (N + j — 2) A. 

Proof. The proof of this lemma is similar to Selberg's original proof for the one-dimensional 
Selberg integral p. The integral is symmetric under relabeling (ni,mx), . . . , {n N} m N ) — )■ 
(t7 p (i), m p (i)), . . . , (n p (jv), m p (jv)) for some permutation of the indexes 1, . . . , iV. The sum, 
Ei=i % + = N{N-1)X, follows because [10 is a homogeneous polynomial. The remaining 
two can be found by letting Xj = yj in [10] and applying the proof for the one-dimensional 
case jfj. Since the proof is only about mj + nj and not xj and yj, those proofs apply here. □ 



Now, consider when we break the 2D Selberg into a sum of products of integrals, 

<Sjv(Ai, A 2 , A3, A4, A) : 
where 

s p ,at(v 1 Ai, v 2 A 2 , A) 



II Yl ) a P,<?(A)s P ,jv(v 1 Ai,v 2 A 2 ,A)tg,jv(v 3 A 3 ,v 4 A4,A), (11) 

a<k<l<N p k i=0/ 



[0,1]' 



X 



n k-^s, 



(12) 



l<fc<«<7V 



and 



t^(v 3 A 3 ,v 4 A 4 ,A) = / flyJ^il-x,) 



t> 4 A 4 



(13) 



l<fc<Z<iV 

where and gjy are symmetric matrices with p^ = = and otherwise pki + to = A, 
n = J2k<lPki> N(N - l)A/2 - n = J2k<ilM> vk = ( u i> ■■■> v i) where the superscript is an 

- 1 = v 2 = „ 
j j j j 

r(i + A) 



index rather than a power, and v) = v 2 = t> 3 = v\ = 1. Also, 



°/mz(A) 



r(i+p w )r(i + 9w ) 



We also have the sum: 



s P ,w(v 1 Ai, v 2 A 2 , A) = 2J c «i,-, 

m,...,njv 

r(i + ujAi + n,)r(i + v 2 \ 2 ) 



n N 



X 



A' 



n 



L T(2 + v}X 1 +v]X 2 +n j ) 



(14) 



but unlike in [TUJ here we know that ^\ rij = 2n. A similar formula holds for t q ^. These 
are essential for the following lemma: 

Lemma 2. For < n\ < n 2 < ■ ■ • < n^, the non-zero c ni)tni) ... jnjVjmjv occur when m, . . . ,tin 
satisfy the inequalities n min (j) < n,- and n max (j) > n^, w/iere n min (j) = \Y,k<i<j 2 Pki/ j] > 
^max(j) = Y.k=i 2 Pjk - \Y.k<i<N+i-j 2 Pki/( N + 1 - jj] and M means the smallest integer 
greater than ■. And a similar lemma involves replacing n with m and p with q. 
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Proof. Consider A^ N \xi, . . . ,x N ) = Yli<k<i<N( x k ~ %i) 2pw , which is a homogeneous poly- 
nomial of degree J2i<k<i<N 2 PM- Based on the ordering, n N > \J2i<k<i<N 2 Pki/N] . 
Since is divisible by for each j = 1,...,N — 1, rij > \Yli<k<i<j^Pki/j1' 

Now, consider that (xf. — x{) 2pkl = xl m x 2pkl (l/xk — l/xi) 2m , hence A^ N \xi, . . . ,x^) = 
lljLi^f i=l * W A^(l/xi, . . . , l/x N ). The allowed powers on the r.h.s., say, < n[ < n' 2 < 
• • • < n' N are given by the powers on the l.h.s. as n'j = J2f =1 2pij — rtM+j-i and must satisfy 

n'i> IEk^W/I- □ 

Lemma 3. The expression[TQ can be evaluated as a function of Xj to give, 



Sp.AT^Ai, V 2 A 2 , A)t 9iA r(v 3 A 3 , V 4 A 4 , A) = C PiN C q!N (X) J | 

r(l + ^ 3 A 3 + m min (j))r(l + ^ 4 A 4 + m min (j)) 
r(2 + y 3 A 3 + <A 4 + m max (j)) 



V r(l + vjXi + n min (j))r(l + v)X 2 + n min (j)) 
._ i V{2 + v)X 1 + v 2 X 2 + n m ^{j)) 



x 



3 ° j 

Proof. From the inequalities for rij, 

r(v}Ai + 1 + ^ 



T(2 + ujAj + v 2 A 2 + n 3 -) 
T(l + v)X 1 + n min {j)) 



T{2 + v)X l + v 2 X 2 + n^{j)Y n ^^ 
where q nj (VjXx,VjX 2 ) is a polynomial of degree n max (j) — % in A 2 . Thus 

v r(l + A 1 + n i )r(l + A 2 ) 



n 



r(2 + Ai + A 2 + n j ; 



3=1 

N 



rtf i\ ^ , 1 - r r(l+t; 7 1 A 1 +n min (j))r(l + ^A 2 ) 
Q (v 1 X 1 ,v z X 2 ) [I 



^ r{2 + vjX 1 + v]X 2 + n max (j)) 
where the polynomial Q'( vl Ai, v 2 A 2 ) has degree 

N 

^2n max (j) - rij. 

3=1 

in A 2 . 

The expression O is a linear combination of terms of this form. Therefore, 

v T(l + v)X 1 +n mhl (j))ni + v 2 X 2 ) 



Sp,N 



7=1 ^ 



+ tj)Ai + v 2 X 2 + (j)) 



gVAi, v 2 A 2 ) -pr r(l + VjXi + n min (j))T(l + v 2 X 2 + n min (j)) 
^(v 2 A 2 ) Al r (2 + t;jA 1 + t; 2 A 2 + r2 max (j)) 



(15) 



g^A^fA,) (16) 
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where Q'(v 1 Ai, v 2 A 2 ) is a polynomial of degree at most ^jLi n maxO) — n j and 

^ (v As) = M r(i + , 2 a 2 ) 

is a polynomial of degree X]j n min(j) which, by Lemma [2] is J2j n ma,x(j) — Uj- Since s Pi at is 
symmetric in v x Ai and v 2 A 2 (by a change of variables Xj — > 1 — sc 3 -), 

g(v% 1 v^g) _ QCv^.v'Ai) 
i?(v 2 A 2 ) i2(v!Ai) ' 

The right hand side is a polynomial in A 2 hence Q(v 1 Ai, v 2 A 2 ) is divisible by _R(v 2 A 2 ), but 
the maximum degree of the two is the same. Therefore, ^t^^v^ must be independent 



of A 2 . By symmetry it is also independent of Ai and can only depend on TV and A 6(. An 
identical proof can be done for t q>N using A 3 and A 4 . Note that, because these two proofs 
can be done separately, the relabeling needed for the initial step can be done first for s P) n 
and next for t q ^ since the integrals are separable. □ 

Lemma 4. Let pj = Yl!k=iPjk an d Qj — Ylfk^Qjk- The limiting case 

lim (1 + Ai)s PjiV (v 1 Ai, v 2 A 2 , A) 

Ai— 1+ 

= Ns p , N ((2p N1 - 1, ... , 2p N , N „ 1 - 1), A 2 , A) (17) 
tfPm 7^ 0. Ifpm, ...,PNj = andp 12 , ...,p lj = 0, 
lim (l + A^+^v^v^A) 

Ai— 1+ 

3 i 
= (N - j)s p>N ((2(p Njj+1 + ^2p j+h i) - 1, . . . ,2(pAr i7V _i + ^2pn-i,i) - 1), A 2 , A) (18) 

i=i i=i 

and likewise for t q ^ replacing p with q. 

Proof. Since the integral is symmetrical in x±, yi, . . . , xn, Vn and 

s p (Xi, X 2 , X) — N\ / dx jy / dxjy-i ■ ■ ■ / dx± x 

Jo J xn J %2 

N 

1=1 j<k 

The stated result follows from the general principle, 

lim (1 + a) f drr a f(r) = /(0). 

a->-l+ J Q 

7 



except when pjvi, • • • , Pnj = and pu = for all k = 1, . . . , j and I = 2, . . . , j + 1, in which 
case the limit blows up, and we use the following formula: 

s p (Xi, A2, A) = (iV — j)\ / (ixAr / dxi ■ ■ ■ / dxj / dxN-i--- / x 

N 

]Jx^(l - x^Hixj - x k ) 2 ^. (20) 

J=l j<k 

The limit is identical for t q replacing p with q. □ 

Combining Lemmas [3] and H] gives the recurrence relations in[6j Note that the formula is 
finite when any of the arguments to the T functions is 0. So, the full formula is, 

A 



n e 

K l<k<l<N p H =0/ 



r(l + A) 

Cp,N c q,N\ 



r(i+p y )r(i + 9o -) 



JL r(i + Ai + E fc <i<j Wig + ^2 + E fc<i<J - Wi) 

11 r(2 + Ai + A 2 + 2 Pi - E*<,sY + i-i 2 p«/(^ + 1 - j)) 

r(i + a 3 + E fc <K, 2Wi)r(i + a 4 + 2Wj) 



X 



(21) 



T(2 + A 3 + A 4 + 2 Qj - Zk<i<N + i-i Zquil \N + l-j)Y 

where q^i = X—pki- This formula, while complex, is computable in polynomial time. The sum 
is iV(iV — l)A/2 terms and there are N 2 products in each term. Therefore, it is computable 
in 0(N 4 ) time. 

Substituting this into the equation in Lemma |3] gives the evaluation 0] of the two- 
dimensional Selberg integral for all non-negative integers A. For = 2 this formula gives 
Ginibre's result [8], which, until now, was the only successful integration of this integral for 
all N. 

Finally, the 2D, one component plasma's partition function (j4]) is easily found using 
Stirling's formula T(l + n) ~ \ / 2nn(n/e) n . 

The main results of this investigation are the exact integration of the 2d one-component 
plasma for even /3, which has been an open problem (in the general form) since its 1965. It 
is likely that this result may apply to a number of other results related to non-Hermitian 
random matrices as well as vortex statistics and plasmas (such as the guiding center plasma) 
such as are mentioned in the introduction. Perhaps most surprising about the 2D Selberg 
integral is that the formula is essentially a sum of products of ID Selberg-style integrals 
(albeit more general). This raises the question as to how many other statistical results 
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applied to one-dimensional ensembles could also be applied to 2d ensembles or even 3d 
ensembles (e.g., quaternion matrices). It may be that a number of such results in interacting 
particle theory could be found using the Selberg integral or other solved Id integrals for large 
dimensional ensembles. 
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